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Abstract 

Four-dimensional simplicial quantum gravity is modified either by coupling 
it to U(X) gauge fields or by introducing a measure weighted by the orders 
of the triangles. Strong coupling expansion and Monte Carlo simulations are 
used. Although the two modifications of the standard pure-gravity model are 
apparently very distinct, they produce strikingly similar results, as far as the 
geometry of random manifolds is concerned. In particular, for an appropriate 
choice of couplings, the branched polymer phase is replaced by a crinkled phase, 
characterized by the susceptibility exponent 7 < and the fractal dimension 
du > 2. The quasi-equivalence between the two models is exploited to get 
further insight into the extended phase diagram of the theory. 



1 Introduction 

As this work follows a couple of other papers which we have devoted to the study of 
the phase diagram of 4d simplicial gravity [|], ^] , we shall not expand much on our 
motivations. Instead, in order to save space, we shall enter without further ado into 
the main body of the work. Let us only observe that following the surprising dis- 
covery (l], H| that the phase transition between the so-called crumpled and branched 
polymer phases in 4d simplicial gravity is discontinuous, several groups have tried 
to figure out what kind of modification could render the theory more realistic. One 
track, explored up to now in 3d only |||, |j| , consists in modifying the measure in the 
partition function following a recipe proposed some time ago in Ref. Q. Loosely 

1 Laboratoire associe au C.N.R.S. 
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speaking, this corresponds to the introduction of a R 2 term in the actionF). Another 
idea consists in introducing matter fields 0. Indeed, in the continuum formalism, 
one can argue that adding conformal matter fields in 4d has the effect similar to 
that obtained by reducing their number in 2d || and might therefore bring one from 
the branched polymer to a physical phase (analogous to the Liouville phase in 2d). 

In Ref. Q the effect of coupling non-compact U(l) gauge fields to 4d simplicial 
quantum gravity was studied. It was found that for more than two gauge fields 
the back-reaction of matter on geometry is strong, that the degeneracy of random 
manifolds into branched polymers does not occur and that there is some evidence 
for a new "smoother" phase of the geometry, to be called hereafter the crinkled 
phase. This apparently confirmed the speculation put forward in ||. The results 
were obtained from the strong coupling series, whose introduction in this context 
has been quite novel, and from Monte Carlo simulations with lattices of relatively 
modest size. An eventual extension of the study to larger systems was announced. 

We have extended our Monte Carlo simulations to systems with up to 32K 
simplexes, for 3 copies of gauge fields. We have also calculated more terms of 
the strong coupling series. While our results are still consistent with the scenario 
described above, we do observe some disturbing inconsistencies in the behavior of 
different geometric observables which, taken at face value, seem incompatible with 
the existence of a single phase transition at a finite value of the Newton coupling 

We have noticed, on the other hand, a curious quasi-equivalence between 
the model with gauge fields and that with a properly modified measure. We have 
concentrated on this issue, achieving a better understanding of the results obtained 
earlier, and we have explored further both the phase structure of the model and the 
nature of the hypothetical crinkled phase. All these developments will be described 
in the following sections. 

2 Gauge fields on a random manifold 

The action is a sum of two parts. The first is the Einstein-Hilbert action, which for 
a 4d simplicial manifold reads: 

S G = -k 2 N 2 + K4A4 , (1) 

where denotes the number of fe-simplexes. The second part is (for one copy of 
the gauge field): 

s M = Y1 °(*«6c) KU + Mhc) + A(l ca )f , (2) 

where A(l a b) is a non-compact U(l) gauge field living on link l a b and A(l a b) = —A(lf, a )- 
The sum extends over all triangles t a b c of the random lattice and o(t a b c ) denotes the 

2 In this text the expression "R 2 term" refers to any combination of terms quadratic in Riemann, 
not only to the square of the scalar curvature. For early simulations with R 2 term in the action see 
also@|. 
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Figure 1: The specific heat, CV, and the susceptibility of the order of the most singular 
vertex, Xp > f° r three copies of gauge fields coupled to simplicial gravity. This is for N4 up 
to 32K. 



order of the triangle t a b c , i.e. the number of simplexes sharing this triangle. We 
introduce / copies of gauge fields. 

We work in a pseudo-canonical ensemble of (spherical) manifolds, with almost 
fixed A^4. The model is defined by the partition function: 

Z(k 2 ,N 4 ) = Y J W(T) f 'j[dA{l) e-Sc-SM-SiNi-Ni) 2 _ (3) 
T ^ leT 

The sum is over all distinct triangulations T and W(T) is the symmetry factor 
taking care of equivalent re-labelings of vertexes. The prime indicates that the zero 
modes of the gauge field are not integrated. As is well known, one must allow the 
volume to fluctuate. The quadratic potential term added to the action ensures, for 
an appropriate choice of 5, that these fluctuations are small. For more details see 
Ref. g. 
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The result of our Monte Carlo simulations of the model defined by Eq. (|3|), 
for 3 copies of gauge fields, is summarized in Figure 1. There we plot both the 
specific heat: Cy = k?,((Nq) — (No) 2 ) /N4, and the susceptibility of the order of the 
most singular vertex pv. Xpo = ((Po) ~ (Po) 2 )/-^4- This is done for N4 up to 32K. 
The specific heat has a peak at K2 ~ 2.5; the location of the peak is stable but its 
height appears to saturate as the volume is increased. This could indicate either a 
continuous phase transition, of 3rd or higher order, or simply a cross-over behavior. 
The susceptibility Xp a l so has a peak, rising with the volume, signaling a change 
in the geometry. Its location, however, moves towards larger and larger values of 
the coupling constant k%. As collecting data at 32K already required a considerable 
effort, going to a significantly larger volume is, for the moment, out of question. 
Thus we cannot tell whether the the location of the peak will tend to a finite value 
of «2 or will go to infinity. 

The apparent inconsistency in the behavior of those two geometric observables 
- their fluctuations seem to be uncorrelated — is rather worrying. A pessimistic 
view is that this implies the absence of a true phase transition in the model; there 
will be no singularity in the infinite volume limit, and the differences in scaling 
behavior of the baby universe distribution and po in comparison to the crumpled 
phase, which we have observed and which indicate a new phase, disappear in the 
thermodynamic limit. While our simulations cannot rule this out, another plausible 
explanation emerges as the phase diagram is explored in more details as we do in 
the next section. We will return to this discussion in the last section. 



3 Modified measure 

The physical variables in the partition function Eq. (|3|) are not the gauge fields, 
but the plaquette values. Let us replace in our model, for the moment without any 
justification, the integration over fields by the integration over plaquettes, which are 
by the same token assumed to fluctuate independently. They can then be integrated 
out, giving rise to a measure factor 

N 2 

M(T) ~ Ho(t^ (4) 
j=i 

with (3 = —f/2. This suggests to compare the two models, the one defined in the 
preceding section and that with the measure factor given by Eq. (||). We have done 
that leaving the parameter [3 free, since it is likely that the estimate (3 = —f/2 is 
correct, if at all, in the limit / — » 00 only. 

First we use the strong coupling series. We have extended the calculation 
presented in Ref. @0 up to iV 4 = 38. The new results are summarized in Table 1. 
In Figure 2 we show a plot of the susceptibility exponent 7 versus f3, calculated in the 
large-K2 phase («2 = 10) using the ratio method [10]. The values of 7 are extracted 



3 There are typos in Table 1 of Ref. ||: The number of triangulations of volume 24 with 9 
vertexes is 34, not 13. Also, the weight Wo of triangulations of volume 28 with 11 vertexes is 
77057 1 = 539400/7, not 77057. 
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Table 1: The number of different graphs N g , for volumes N 4 = 32, 34, 36 and 38, and the 
corresponding weights Wf(N±, No). This is shown both for pure gravity, and one and three 
gauge fields coupled to gravity (/ = 0, 1, and 3). All weights are normalized with the value 
at N 4 = 6. 
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assuming a large-volume behavior of the canonical partition function: Z{N^) ~ 
exp^A^A^ -3 , where \x c is the critical cosmological constant. We observe a range 
of (negative) couplings /3 where a reliable estimate of 7 is obtained, converging as 
more terms of the series are included in the analysis. And, just as for the gauge 
field model, in this interval of (3 the exponent 7 decreases continuously from the 
branched polymer value, 7 = 1/2, to a large negative value. 

The similarity of results obtained with the two models, gauge fields and mod- 
ified measure, is in fact more spectacular. In Figure 2 we also trace 7 versus /, 
making the identification 

" = -14 < 5 > 

In an interval of the values of f3 the two curves coincide! Actually, there is close 
agreement in the parameter range where the ratio method seems to give reliable 
values of 7. Note, that the leading term in Eq. (|5|) coincides with our earlier estimate. 
Of course, Eq. (|5|) is meaningless for too small values of / (— /3): there 7 = 1/2 for 
both models and the mapping is trivial. 

We have further checked this remarkable equivalence in a numerical exper- 
iment, measuring the effective change to the action Eq. (|l|) stemming from these 
two seemingly unrelated modifications of the standard model. In a Monte Carlo 
simulation of pure gravity, with K2 = f3 = f = and N4 = IK, we measured the 
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Figure 2: Variations of 7 with both the coupling (3 to the measure term Eq. (||) (solid line), 
and with the number of gauge fields / (dashed line), calculated from the series expansion. 
The value of / has been rescaled using the relation Eq. (||). Also shown are values of 7 
measured in Monte Carlo simulations with the modified measure for K2 = 3, varying 0, on 
volume N4 = AK . (Insert: The correlations between the effective actions stemming from 
integrating out the gauge fields and from the measure term, Om and Ot-) 

effective actions: 

M = JjlogCofo)) and O v = - log I J, (6) 

3 ' 

where A is the determinant associated with the integration over one species of gauge 
fields. We observe a very strong linear correlation between those two observables, as 
shown in the insert in Figure 2. A linear fit yields a slope —0.65, roughly compatible 
with Eq. (|). 

Further corroboration comes from Monte Carlo simulations at j3 = —1.75 
(corresponding to / = 3 in Eq. (||)), for N± = 4K and 8K and for K2 ranging from 
to 4.5 . We do not wish to drown the reader in figures. It suffices to say that, up 
to a shift in ki not larger than ~ 0.5, the results concerning the first two cumulants 
of Nq and the orders of the three most singular vertexes are practically the sameQ. 

As simulating the model with a modified measure is less CPU-demanding than 
working with several copies of gauge fields, it has been possible to explore the phase 
diagram in more details. In addition, the coupling (5 can be varied continuously, 
whereas we are restricted to integer values of /. We did simulations for K2 = 3, 
JV4 = 4K, and for several values of 0. We observe strong fluctuations in geometry 
both at (3 w —1.5 and /3 ~ — 5. This is indicated by peaks in the susceptibility of 

4 The shift in K2 depends, however, on the observable. 
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Figure 3: The fluctuations in the order of the three most singular vertexes, po, p% and P2, 
for simplicial gravity with a modified measure. This is for K2 — 3 and volume is N& = AK . 

the various geometric observables; in Figure 4 we show this for the orders of the 
three most singular vertexes. The first signal corresponds roughly to the transition 
from branched polymers to the crinkled phase, as indicated by the series expansion, 
whereas the latter is a transition to a crumpled phase for (3 > — 5. At the moment 
we have only investigated this at one lattice volume, obviously further exploration 
is needed to establish that this indeed corresponds to two distinct transitions. 

To end this section let us mention that this universality in the back-reaction 
on the geometry also holds for other modifications of the standard model. Using 
the series expansion we have investigated the effects of adding: (a) Gaussian scalar 
fields, (b) a discretized R 2 -term as used in Ref. Q, and (c) a modified measure 
using a product of vertex orders. In all three cases it is possible, by an appropriate 
rescaling of the corresponding couplings, to map the extracted values of 7 on the 
curve in Figure 2, although the agreement is not as spectacular as that between the 
model with gauge fields and the one with a measure modified according to Eq. @), 
respectively. 

4 The nature of crinkled manifolds 

We have explored further the nature of the hypothetical crinkled phase in Monte 
Carlo simulations. We have calculated the exponent 7, from the distribution of baby 
universes, along two lines in the phase diagram: l\ = {K2 = 3, (3} and I2 = {«2,/3 = 
—3.5}. The values of 7 measured along l\, included in Figure 2, agree reasonably 
with the predictions of the series expansion. The values of 7 measured along I2 are 
shown in Table 2. For small values of K2 it is not possible to extract any reliable 
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Table 2: Values of the string susceptibility exponent 7, measured in Monte Carlo simula- 
tions at N4 = 4K , for f3 = —3.5 and varying k 2 . The series expansion predicts 7 ~ —6.2, 
for /3 — —3.5 and large values of K2- 
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Table 3: The extracted values of the spectral dimension d s for simplical gravity with a 
modified measure. This is for K2 = 4.5 and both for (3 = (branched polymer phase) and 
for j3 = —1.75 and -3.5 (crinkled phase). 
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value. This is to be expected as the model is in the crumpled phase. For k 2 > 3.5, 
however, we get a consistent value 7 ~ —4.5, compared with the prediction of the 
strong coupling expansion, 7 —6.2, obtained for f3 = —3.5 and n 2 large. 

Other exponents characterizing the fractal geometry are the HausdorfF and 
spectral dimensions of the manifolds, dn and d s . The former is related to the 
volume of space within a sphere of geodesic radius r from a marked point: v{r) ~ 
r dn , whereas the spectral dimension defines the return probability for a random 
walker on the triangulation: p(t) ~ t~ ds / 2 [|ll|. The time t is measured in units of 
jumps between neighboring vertexes, with hopping probability given by the inverse 
of the coordination number^]. On smooth regular manifolds those two definitions 
of dimensionality coincide. However, on highly fractal manifolds, like the ones that 
dominate the partition function Eq. (||), they are in general different. 

We have extracted the Hausdorff dimension for K2 = 4.5 and P = —1.75 and 
—3.5, from the expected scaling behavior of the average distance between two sim- 
plexes: (rij)N 4 ~ N\ H . Measurements at volume ./V4 = 4K to 32K were included 
in the fit. For (3 = —1.75 we got dn = 3.57(16), which should be compared to 
d H = 3.97(15) quoted in Ref. |] for / = 3 and k 2 = 4.5. For = -3.5 we got 
djj ~ 5, but this estimate is less reliable. 

The measured values of the spectral dimension at k 2 = 4.5 for ./V4 ranging 

In calculating both the Hausdorff and spectral dimensions we used the dual graph. This is 
more natural as we measure at fixed volume N&, i.e. fixed number of vertexes in the dual graph. 
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from 2K to 16K and for three values of (5 are shown in Table 3. In the branched 
polymer phase, at (3 = 0, we get d s « 1.33, to be compared to the theoretical value 
for generic branched polymers: d s = 4/3 [O]. In the crinkled phase, on the other 
hand, we get values for the spectral dimension significantly larger than 4/3 and 
which, moreover, seem to increase as (5 is decreased: d s ~ 1.5 at f3 = —1.75, and 
d s ~ 1.75 at (3 = —3.5. In all cases, the values obtained at different volumes agree 
within the numerical accuracy. 



5 Discussion 



The quasi-equivalence between the model with gauge fields and that with modified 
measure cannot be a coincidence. The simplest explanation is that the correlation 
between plaquettes falls rapidly with the distance on the lattice. This sheds a 
new light on the results of Ref. |Q. The effective action used in the speculation of 
Ref. j|] is derived from trace anomalies, assuming conformal invar iance. This regime 
is apparently different from the one we observe on our disordered lattice. If so, then 
the mechanism of suppression of polymerization observed in the model is not the 
one which has been expected and a faithful implementation of the idea of Ref. |^] 
remains an open problem. 

We are now in a better position to discuss certain points which were left 
obscure in Ref. Q. In particular, mean-field arguments give perhaps some insight 
into our results: 




Following Ref. 12] define a mean-field model by the canonical ensemble par- 
tition function 



5N 4 , a < -2, (7) 



where P denotes partitions of the integer 5N4 into No "vertex orders" q and a is an 
adjustable parameter^. At large N± one can write 

z ~ [*dre N *l Kr + fW , r = N /N 4 , (8) 
Jo 

the function f(r) having two remarkable properties: f'(r) =const (call it —K cr ) for 
r < r*, and f"(r) < for r > r*. Furthermore r* increases with —a (see Ref. |l2j 
for details). 

Thus, in the thermodynamic limit and as long as r* < 1/4, the model has 
the following two phases: for k < K cr one has r = 0, while for k > K cr one finds 
r = r sp, where r sp > r* is the position of the saddle-point. The "latent heat" at the 
transition equals r* . Loosely speaking, the two phases correspond to crumpled and 

6 It might appear more natural here to assume that triangle and not vertex orders are the 
variables fluctuating independently (up to the global kinematic constraint). It turns out that one 
would then obtain an unphysical model, predicting singular triangles in the crumpled phase, in 
variance with observations. 
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branched polymer geometries, respectively: For k < K cr a few vertexes have orders 
~ A4, while for k > n cr vertex orders are bounded. 

The situation changes when r* > 1/4: there is no saddle point in the integra- 
tion range and at k = K cr the system jumps from r = to r = 1/4 (this argument, 
implicit in Ref. fl~2[| , is emphasized in Refs. [b|]). In the model, there are still 
singular vertexes (i.e. with order ~ N4) when k > K cr , but the order of the most 
singular vertex drops suddenly as one moves across k = K cr . Hence, at the 
corresponding susceptibility is expected to have a peak with height oc N4. 

The overall picture shows some similarity to what we observe, especially if one 
recalls the results mentioned at the end of Section 3, implying that the increase of 
—a can mimic the increase of / or —/3 in the modified simplicial gravity model. 

Although the mean field model clearly helps understanding our results, it 
should perhaps not be taken too literally. In particular, the model predicts that the 
only coupling space region where the crinkled phase survives in the thermodynamic 
limit is the region where Nq/N^ = 1/4, i.e. where the naive Regge curvature sticks 
to its upper kinematic bound. This scenario is plausible. If it is true, the crinkled 
phase is an unlikely candidate for a physical phase of quantum gravity. In the data 
(see Figure 1) the transition point seems to run towards large values of K2, where 
indeed {Nq)/N^ — > 1/4. However, one observes Nq/N^ = 1/4 also in the branched 
polymer phase, at finite volume and large k,2- This could be just a finite size effect. 
We have no evidence for a dramatic jump of (Nq)/N4 : as the system moves from the 
crumpled to the crinkled phase. Actually, the behavior of the specific heat shown 
in the upper part of Figure 1, seems to exclude any 1st order transition. Also, the 
most singular vertex susceptibility increases slower than the volume. 

We should also mention another point where the mean field model appears in 
variance with data. As already mentioned, the qualitative prediction of the model 
is that the latent heat increases with —a (or — (3, or /). However, in Ref. [|| we have 
noted that the latent heat at the transition point decreases (by a factor of 2 at 32K) 
as one moves / from to 1. This observation is compatible with the claim made 
in Ref. ||, but in 3d, that the transition becomes of second order when the power 
in the measure (in their case the power of the vertex order) becomes sufficiently 
negative. 

We believe that the results of our study can be tentatively summarized by the 
phase diagram^ of Figure 5. There are three distinct regions corresponding to the 
branched polymer phase, the crumpled phase and the hypothetical crinkled phase, 
respectively. The solid line BA represents the line of phase transitions separating 
the crumpled and the branched polymer phases. The two dashed lines, BC and BD, 
separate the crinkled phase from, respectively, the crumpled and branched polymer 
ones. It is unclear yet whether these lines represent genuine phase transitions, in 

7 This phase diagram suggests an explanation for the inconsistency observed in Section 2 in 
the simulations with 3 gauge fields. As discussed in Section 3, the transition from the branched 
polymer phase to the crinkled phase occurs at /3 ~ —1.5. This is accompanied by strongfluctuations 
in geometric observables such as the maximal order of vertexes. By the relation Eq. (pi) this value 
of /3 corresponds to / ~ 2.5. Hence simulating 3 gauge fields might have been an unlucky choice. 
It it possible that varying K2 we have followed the line BD in the phase diagram, uncomfortably 
close to the transition region. 
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Figure 4: A schematic phase diagram of simplical gravity modified with the measure Eq. 
At least along a portion of the line BA the phase transition is discontinuous. The same 
phase diagram should hold for gauge fields coupled to gravity, Eq. (||), replacing (3 by the 
number of fields / rescaled by the relation (||) . 



which case the point B is a tricritical point, or a cross-over behavior. Further 
simulations are needed in order to clarify this issue and to verify that this phase 
structure survives in the thermodynamic limit. 
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